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ABSTRACT 

Characteristics of the initial condensate in the recent experiment on Bose-Einstein con- 
densation (BEC) of ®^Rb atoms in an anisotropic magnetic trap is discussed. Given the 
aspect ratio R, the quahty of BEC is estimated. A simple analytical Ansatz for the initial 
condensate wave function is proposed as a function of the aspect ratio which, in contrast 
to the Baym-Pethick trial wave function, can be used for any interaction strength, repro- 
duces both the weak and the strong interaction limits, and which is in better agreement 
with numerical results than the latter. 
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1 Introduction 



Bose-Einstein condensation (BEC) is a phenomenon where a macroscopic number of par- 
ticles is in the ground state of the system at finite temperature. The phenomenon of 
BEC plays significant roles in many branches of physics Because of the presence of 
strong interactions, BEC has been inferred rather than directly observed so far. Recently, 
however, three different groups have reported the direct evidence of BEC in weakly in- 
teracting systems of atoms such as rubidium p|, lithium [Q, and sodium f^, confined in 
anisotropic magnetic traps and cooled down to very low temperatures. These experiments 
show promise of becoming a new laboratory for quantum statistical phenomena that are 
inaccessible to other conventional techniques and of enabling us to experimental study 
phenomena that have been addressed only theoretically, such as spontaneous symmetry 
breaking and decay of unstable macroscopic states. They may also advance our under- 
standing of superconductivity and superfluidity in more complex systems. Moreover, the 
technology used in the experiments has the possibility to be extended to produce a verita- 
ble atomic laser that is bound to have many applications in pure science and technology 

In recent experiments in the system of rubidium atoms and sodium atoms the 
onset of BEC is signalled by a narrow peak on top of a broad thermal velocity distribution 
centered at zero velocity. This peak exhibits the nonthermal, anisotropic velocity distri- 
bution expected of the minimum-energy quantum state of the magnetic trap in contrast 
to the thermal, isotropic velocity distribution observed in the broad uncondensed fraction. 
The parameter which characterizes the asymmetry of the velocity distribution function 
is the so called aspect ratio R = \J {pD / {pi). In the experiment of Anderson et al. 0, 
rubidium atoms are initially trapped and cooled down in a strong magnetic trap which 
can be described as a three-dimensional (3D) harmonic potential cylindrically symmetric 
about the 2;-axis, with tunable frequency uo^ (in the z direction) and ujx_ = Uz/X (in the 
xy-plane), with the asymmetry parameter A = -\/8- The corresponding oscillator lengths 
are a_L(^) = (/i/mcj^(^))i/2 ^ 1.25(0.74) x 10" '^cm, where m is the atomic mass. After 
some time, the cloud of rubidium atoms is adiabatically released to a weaker magnetic 
trap whose spring constants are 10 times weaker than when BEC forms. The condensate 
is then examined after ballistic expansion from the weak trap. The ballistic expansion 
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is properly modelled numerically by a self-consistent wavefunction calculated by Holland 
and Cooper H which generalizes previous investigation of the symmetric evolution M. 



Similarly to previous recent studies |T0[, we shall confine ourselves to the system 



of ^^Rb atoms and we will discuss the characteristics of the initial condensate formed. 
Although the characteristics of the initial condensate have not been measured directly 
yet, there is hope that it will be possible so in future experiments. 

A theoretical picture of the initial condensate was produced within the Hartree-Fock 
(HF) approximation p, ^ |T^ using the Ginzburg-Pitaevskii-Gross (GPG) energy func- 
tional and associated with it the Nonlinear Schrodinger Equation (NSE). Excitations, 
using the technique of the Bogoliubov transformation, have been described by Fetter . 
Baym and Pethick IQ have gained an insight into the problem by assuming that, similarly 
to the noninteracting case, a gaussian gives a reasonable variational ground-state wave 
function, the only effect of interactions being a renormalization of oscillator frequencies. 
They show that the first effect of interactions is to reduce the density of the cloud of par- 
ticles in the central region from the free particle situation and expand it in the transverse 
direction. These qualitative features were confirmed by Edwards et al and Dalfovo and 
Stringari |0 by solving the NSE numerically. 

In the present paper it is shown that although, qualitatively, the Baym and Pethick 
(BP) scenario [|| is correct, nevertheless, the BP wave function does not describe BEG 
regime well. As we show in the next section, this can be explained by the fact that the 
BP wave function is actually the square root of the first order density function in the 
high temperature expansion of the partition function of the 5-function interacting Bose 
gas in the kinetic energy, and hence describes high temperature properties of the system. 
Higher order corrections are needed to obtain low temperature properties such as BEG in 
agreement with the numerics. In particular, this explains the unreasonably high aspect 
ratio in BP estimations (up to 4.2, i.e., 2.5 higher than in the noninteracting case). All 
these prompt a search for another trial variational function for the ground state. In the 
section |^ we derive a simple analytical Ansatz [see Eq. (|^) and Fig. |I[ which, for a 
given numerical value of the aspect ratio, describes well the ground state properties of 
the system for all values of the interaction strength. The Ansatz interpolates smoothly 
between the weak and the strong interaction case. It is worthwhile to notice that the 
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Figure 1: A comparison of the x-dependence of the numerical solution of NSE for the 
ground state [10] (dashed line) and our approximate solution (solid line) in the case of 
= 5000 atoms of ^^Rb, when A = 502 and Cansatz = 2.2. The value of the aspect ratio, 
R{A, A) = 2.3, is taken from Ref [10]. 

latter case cannot be decribed by the BP wave function and, instead, it is described by 
the Thomas- Fermi approximation in Ref. [^. Using this Ansatz, correlation effects can 
be considered |]T^ . 



2 Baym-Pethick trial wave function and high tem- 
perature expansion 

Baym-Pethick argued [§] that the initial condensate wave function in the region of weak 
and up to intermediate interactions can be well approximated by a gaussian. Let us take 
a different point of view and ask under which condition this gaussian-like profile of the 
initial condensate wave function can be actually derived. 

If interactions between atoms are neglected, the physical system is equivalent to the 
system of noninteracting harmonic oscillators. In the latter case, one can show that the 
aspect ratio in thermal equilibrium is a monotonically increasing function of the inverse 
temperature /3 = 1/T, 



R{X,P) 



(1/2) + (e/^^"^ - l)-i 

\ (1/2) + (6/3-^ - 1)-1 ' 
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At high temperatures i? ~ 1, and for low temperatures R ~ V^- In the latter case, the 
dominant contribution to the aspect ratio is given by the ground state of the system and 
reflects its anisotropy. 

In the presence of the interactions, the Hamiltonian of the system can be written as 



H 



2m 



dV 



(1) 



where = + y^, a± and are oscillator lengths, and / is the s-wave scattering 



length p, |1^. Qualitatively, the temperature dependence of R{X, f3) preserves the main 



features of the noninteracting case. At high temperatures, the interaction is irrelevant 
and R ^ 1. At low temperatures (for sufficiently small fraction of atoms out of BEC), 
the HF approximation is justified and the ground state wave function (giving the main 



contribution to the aspect ratio) satisfies the NSE. After rescaling of variables JTO|, the 
NSE can be written as 



(2) 



Here, A = SttIN/ a_L characterizes the interaction strength, is the number of particles in 



the condensate {A ~ 520 for A^ ~ 5000 [0), and C = fi/ {huj±) > 0, being the chemical 
potential. In the case of large condensate fraction (strongly interacting case, A ^ 1), the 
kinetic term can be neglected M, |10[ and the ground state (normalized to unity) wave 



function is given by the Thomas-Fermi approximation. 
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(3) 



where 2C = [15AA/(87r)]^/^, and G is the Heaviside step function. The aspect ratio, 
R{A, A), is increasing function of A, and the ground state solution (|^) takes on its maximal 
possible value {R = A) among all ground state solutions to the NSE. In the present case 
(A = 78), this means that the maximal effect of interactions is to raise the value of 
the aspect ratio on 67% with respect to the noninteracting case. Moreover, as shown by 
Dalfovo and Stringari |T0[ the aspect ratio for A = 520 (corresponding to A^ = 5000 atoms 
in BEC) is -R = 2.3, i.e., 37% higher than in the noninteracting case. 

Let us now return to the BP wave function In order to simplify the derivation 
of the BP wave function from the high-temperature expansion, let us for a while, such 
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as in Ref. neglect the anisotropy of the oscillator potential |16|. Because the kinetic 
energy of particles in the system is approximately 200 times smaller than the characteristic 
interaction energy |p, it is reasonable to consider the expansion of the partition function 
of the system, /x), in powers of the kinetic term, 

/i) = / r)D^(x, r)e^° ^ / / ^^^^^ = ^ 

•' n=0 x"^ ) n=0 

(4) 



Here 5*0 is the "unperturbed" action. 



5o 



1 dr 2m 



( 2m \ 



with fields satisfying the periodic boundary conditions, (x, /3) = ^/'"'"(x, 0) and "i/'l^^; 1^) = 
ip{x, 0). Expansion (Q) is the high temperature expansion and we are not exactly in the 
BEC regime. However, the HF approximation is avoided. 

In the first term (corresponding to n = 0) of the sum in Eq. (|^) one finds a product 
of partition functions of the anharmonic oscillators {H^}, labelled by the space point x. 
Using the lattice approximation x = /(mi,m2,m3) = Im with the scattering length / 
(the smallest length in the system) being the lattice spacing, one has for field operators 
iIj~^{x) = "^Am/v^ and tplx) = ipm/Vl^, and 



Here, Z"^{j3,fi) is the partition function of the one-site Hamiltonian H„ 

a", 



Hr, 



2mP 

After redefining parameters 



^ / 2 2 2\ 2Tnl 

— (mi + m2 + mg) - -^/x - An 



^^{x)ip{x) + An (jjj^{x)ijj{x)^ 



2ml' 



■ft, P^P = 



2mP 



P, b„ 



— {ml + ml + mj) 



one has 



Z'"(/?, fi) = Yl exp (-47r/3fc2 + + Anpk - pkb„ 



k=0 



Let us assume that the value of the chemical potential (3 fl is negative and order of 10 
(self-consistency of this assumption will be shown below). By substituting T ~ 10~''K for 
the temperature, one has An (3 ~ 12n x 10^, and 



Z"^(/?,/i) = 1 + e 



I3il-f5br, 



l + 0(e 



-2.57rl0"' 



1 + e 



I3il-f3b„ 



(5) 



The resulting partition function, 

m 

leads to the distribution function {iIj~^{x)iIj{x)) of the Fermi-Dirac type, 

1 



1 + e-l3p.+i3b 

m 

with the chemical potential jl to be determined from the normalization condition, 

7V=/ '"^ . (6) 

Let ~ 5000 be the number of particles in the system. Using that l^/a\_ ^ 2 x 10~^°, 
and, neglecting 1 in the denominator of Eq. (^, one gets e^'^ ^ 1.3 x lO^^'^A^/vr'^/^ which 
implies (3jl ~ In (10^^) ~ —16.1, in full accord with our assumption. 

We cannot justify our treatment for low temperatures. Nevertheless, under the as- 
sumption that our expansion holds up to low temperatures, the profile of the square root 
of the density function iIjq^x) of the system turns out to be the BP trial wave function [Q, 

ry.2 \ 



tjjo{x) = {i/j+{x)iIj{x)) = const x exp 



2a 



7.2 ' 



where a±_ = a^/ (^\fp^. The very fact that the BP wave function is the first order result 
suggest that it may not describe BEC well. Therefore, the aspect ratio obtained from 
the BP wave function may not be reliable. Moreover, the profile of the BP function 
differs significantly from the exact ground state wave function calculated numericaly in 



Refs.g |T0|. 

3 New analytical Ansatz 

We shall show that the ground-state wave function /(r) of the system can be well described 
by a simple analytical Ansatz (cf. Fig. |^), 

f = }-W [Aexp(4C - p2 - R\A, \)z^)\ , (7) 

where W{x), defined as the principal branch (regular at the origin) solution to the Eq. 
We^ = X |T3, is a standard MAPLE function |TB|. Constant C in Eq. (0) is to 
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Figure 2: A comparison of the x-dependence of the ground-state solution of Baym and 
Pethick (sohd hue), strong hmit (dot-dashed hne) and our approximate solution (dashed 
line). The same values of the parameters were used as for Fig. 0. 



be determined from the normalization condition, / p{r)dp'x = 1. To approximate the 
ground-state solution to the NSE, the value of R should be supplied from the numerical 
solution |TD[. On the other hand, given experimental value of R, our Ansatz can serve to 
reproduce the profile of the ground state. 

Eq. d^) is a nonlinear equation and, in the present case, no exact solutions are known 
except for the two limiting cases, A = and A/ ^ Af^. In what follows, we shall 
construct our Ansatz to reproduce correctly the two limiting cases and to interpolate 
smoothly between them as interaction changes. Note that in the strong interaction limit 
the Baym-Pethick trial wave-function cannot be used at all and the Thomas-Fermi ap- 
proximation was used in Ref. ^ in this limit. The main point in our derivation is to use 
instead of the 2nd order differential equation (0) a set of first order differential equations 
which reproduce correctly both the noninteracting limit {A = 0) and the strongly inter- 
acting limit (A/ ^ Af^). In this sense, our Ansatz will be exact in the first derivatives. 
In the noninteracting limit, Eq. (0) reduces to the stationary Schrodinger equation for an 
anisotropic oscillator and the (normalized to unity) ground state wave function is 



/(r) 



TT' 



^3/4 



exp 



- ( 



2 I 2 I \ 2 

X + y + Xz 



(8) 
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Figure 3: A comparison of the ^-dependence of the ground-state solution of Baym and 
Pethick (dot-dashed hne), strong hmit (sohd hne) and our approximate solution (dashed 
line) . 

with C = 2 + X and, in agreement with our previous discussion, R = \/A- One notice that 
/(r) satisfies the set of first order differential equations, 

dj=-xf, d2f = -yf, X-^d,f = R-\X)d,f = -zf. (9) 

On the other hand, if A/ ^ Ap-, Af^/X, the (normalized to unity) ground-state wave 
function is given by the Thomas- Fermi approximation [see Eq. (^], and 

AfdJ ^ -xf, Afd^f ^ -yf, X''Afd,f = R-\X)Afd,f ^ -zf. (10) 

Now, first order differential equations (|TU]) can be combined into the variational principle, 
E[f] = lEjP.Pjdh, where Pi = [{1 + Af)d^ + x] f , = [{1 + Ap)d2 + y] f , P3 = 
[R^'^{A, A)(l -|- Af'^)d3 + z] f. The actual form of P3 in the asymmetric case is fixed by 
the requirement to make simultaneous integration of the first order differential equations 
([TT|) possible. The variational principle implies the following set of first order equations 
for the ground state wave function. 



By integrating Eqs. ( [Tl| ) one obtains that 

/e^/V2 = exp \2C - pV2 - R\A, X)z^ 12 
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from which our Ansatz (|^ follows immediately. Provided that + B?{A^ \)z'^ < AC, f 
can be found explicitly using successive iterations, 

where s = exp[2C — — R'^{A, X)z'^]. Obviously, /(r) given by Eq. (J^) reproduces 
correctly the ground state wave function both in the noninter acting limit {A = 0) and 
in the strongly interacting limit {Af'^ ^ 1), and interpolates smoothly between the two 
limiting cases in the intermediate region (see Figs. |^ and|]). One can verify that if Af'^ ^ 1 

then A / [l/{Apy]{x'^ +y'^ + R^{A, X)z^ + 3Ap)f, and the kinetic term is suppressed 

by the factor (Ap)'"^ with respect to the remaining terms in (|^). 

Our Ansatz can substitute for the BP trial wave function and can play the role of a 
new trial wave function in various variational calculations. Given experimental values of 
the aspect ratio, our Ansatz can be effectively applied to describe the initial BEC wave 
function and to calculate all relevant properties of the initial BEC. We believe that the 
derivation of our Ansatz can be extended to deal with excited states too. 

4 Discussion and conclusions 

In the paper, only the properties of the initial condensate were considered. To find the 
connection with experiment, it is necessary to discuss properties of the system during the 
transition from the strong trap to the weak trap and its subsequent ballistic expansion 
from the weak trap. Obviously, characteristics of the system will change after the ex- 
pansion and will strongly depend on the condition of the expansion (i.e., whether it is 
adiabatic or abrupt). Nevertheless, the present study allows us to give an upper bound for 
the aspect ratio of the condensate after its expansion directly from the strong trap, i. e., 
in the absence of the intermediate weak trap, as it took place in 0]: the aspect ratio of the 
final system is always lower than that calculated for the initial condensate. Indeed, after 
the expansion, (i) there is no more anisotropic potential applied, (ii) the self-interaction 
of bosons, which gives rise to the increase of the aspect ratio with respect to the nonin- 
teracting case, is decreasing. Note that if one can measure the aspect ratio directly after 
the expansion from the strong trap as it should be done in future experiments, it would 
be possible to estimate the real number of the particles in the condensate. 
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We want to emphasize that, as it can be found from comparison of the resuhs ob- 
tained from the BP wave function, the Thomas- Fermi approximation, and from the exact 
numerical solution, the aspect ratio may be a very sensitive characteristic of the wave 
function and may change considerably even when other characteristics are not changed 
by a perturbation (such as an external potential or an interaction). That is why we expect 
that taking into account correlation effects may lead to the considerable changes in the 
aspect ratio although for other characteristics the Hartree-Fock approximation will give 
correct results. We will consider this question in details in a forthcoming paper |T3 . 



Summarizing, given the value of the aspect ratio, both the profile of the ground state 
and the quality of BEC can be estimated. This allows one to estimate the number of 
particles in the initial condensate. We showed that the Baym-Pethick trial wave-function 
(i) is only the first order approximation in the high temperature expansion for the system 
and (ii) does not describe the condensate wave function accurately even for weak and 
intermediate interactions. Note that in the strong interaction limit the Baym-Pethick 
trial wave-function cannot be used at all and instead the Thomas-Fermi approximation 
was used in Ref. In order to describe the ground state of the initial condensate in 
the whole range of the interparticle interactions, we proposed a simple analytical Ansatz 
which, in contrast to the BP trial wave function, reproduces correctly both the weak and 
the strong interaction limit and interpolates smoothly between the two limiting cases as 
interaction changes. 
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Grant GR/J35214, and by the Grant Agency of the Czech Republic under Project No. 
202/93/0689. 

References 

[1] A. Griffin, D.W. Snoke, and S. Stringari, Bose-Einstein Condensation (Cambridge 
University Press, Cambridge, 1995). 



10 



[2 

[3 

[4 

[5] 
[6] 
[7 



[9 

[lo; 
[11 

[12 

[is: 

[14 



M.H. Anderson, J.R. Ensher, M.R. Matthews, CE. Wieman, and E.A. Cornell, Sci- 
ence 269, 198 (1995). 

C.C. Bradley, C.A. Sackett, J.J. ToUett, and R.G. Hulet, Phys. Rev. Lett. 75, 1687 
(1995). 

K.B. Davis, M.-O. Mewes, M.R. Andrews, N.J. van Druten, D.S. Durfee, D.M. Kurn, 
and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995). 

K. Burnett, Science 269, 182 (1995). 

M. J. Holland and J. Cooper, to appear in Phys. Rev. A. 

P. A. Ruprecht, M. J. Holland, K. Burnett, and M. Edwards, Phys. Rev. A 51, 4701 
(1995). 

G. Baym and C. Pethick, Phys. Rev. Lett. 76, 6 (1996) ( |cond-mat/950804q) . 

M. Edwards, R. J. Dodd, C. W. Clark, P. A. Ruprecht, and K. Burnett, to appear 
in Phys. Rev. A. 

F. Dalfovo and S. Stringari, Bosons in anisotropic traps: ground state and vortices, 
to appear in Phys. Rev. A ( |cond-mat / 9 5 1 1 4^ ) . 

V.L. Ginburg and L.P. Pitaevskii, Sov. Phys.-JETP 7, 858 (1958); E.P. Gross, J. 
Math. Phys. 4, 195 (1963). 

A.L. Fetter, Ground State and excited States of a Confined Condensed Bose Gas, 



preprint |cond-mat/9510037 . 



J.M.F. Gunn, K.N. Ilinski, and E. Klepfish, work in progress. 



The s-wave scattering length can be estimated as / ~ lOOoo [0, where oq is the Bohr 
radius (i.e., Z ~ 5 ■ 10~^cm). 



[15] Here the value ~ 5000 is chosen in order to make a comparison our Ansatz with 
the numerical solution of Ref. |]TD[) possible. 

11 



[16] If A 7^ 1, one can rescale coordinate z in Eq. (|ID according \,o z ^ z = \f\z which 
imphes that —d^ + y?z^ ^ A (— 9| + . Therefore, the only effect of the asymmetry 
is that the s-wave scattering length / becomes effectively // A in the z-direction. 

[17] F. N. Fritsch, R. E. Schafer, and W. P. Crowley, Commun. of the ACM, Feb. 1973, 
Vol. 16, No. 2. 

[18] B.W. Char, K.O. Geddes, G.H. Gonnet, B.L. Leong, M.B. Monagan, S.M. Watt, 
Maple V Language Reference Manual ( Springer- Verlag, New York, 1991). 

[19] J.R. Gardner, R.A. Cline, J.D. Miller, D.J. Heinzen, H.M.J.M. Boesten, and B.J. 
Verhaar, Phys. Rev. Lett. 74, 3764 (1995). 



12 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/9512155v2 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/9512155v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/9512155v2 



